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Abstract 

Compound resonances in nucleon-nucleus scattering are related to the discrete spectrum of the 
target. Such resonances can be studied in a unified and general framework by a scattering model 
that uses sturmian expansions of postulated multichannel interactions between the colliding nuclei. 
Associated with such expanded multichannel interactions are algebraic multichannel scattering 
matrices. The matrix structure of the inherent Green functions not only facilitates extraction 
of the sub-threshold (compound nucleus) bound state spin-parity values and energies but also 
readily gives the energies and widths of resonances in the scattering regime. We exploited also the 
ability of the sturmian-expansion method to deal with non-local interactions to take into account 
the strong non-local effects introduced by the Pauli principle. As an example, we have used the 
collective model (to second order) to define a multichannel potential matrix for low energy neutron- 
12 C scattering allowing coupling between the 0^ (ground), 2f (4.4389 MeV), and 0^ (7.64 MeV) 
states. The algebraic S matrix for this system has been evaluated and the sub-threshold bound 
states as well as cross sections and polarizations as functions of energy are predicted. The results 
are reflected in the actual measured data, and are shown to be consistent with expectations as may 
be based upon a shell model description of the target and of the compound nucleus. 
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I. INTRODUCTION 



Predicting the scattering of low energy (E < 10 MeV) protons and neutrons from nuclei 
is of wide-spread interest, as is their capture by same. Intrinsically, the probabilities of 
reactions (measured cross-section data) reflect the inherent structure and modes of excitation 
of the struck nucleus. However, despite the large amount of data so far accrued, there 
remains much that has not been measured, or cannot be, but the knowledge of which is 
prime input to many fields of study; fields as diverse as nuclear-astrophysics/cosmology and 
nuclear stockpile stewardship. 

The most commonly studied process (experimentally and theoretically) is that of elastic 
scattering; the outcome of such data analysis generally being required to begin studies 
of other reaction processes. Typically, low-energy elastic scattering data show resonances 
upon a smoothly varying background, with cross-section magnitudes of the order of barns. 
Resonances can be quite varied in their character. Particularly their widths (full widths 
at half maximum, FWHM) vary from a few eV to over an MeV. For light mass targets 
the resonances in the low energy regime tend to be distinct and but a few in number. As 
mass increases, however, the numbers of resonances rapidly increase and the first tends 
lower in energy. Many strongly overlap, so that a statistical approach to the analysis of 
scattering becomes feasible and utilitarian. But in the main, data analysis to date has been 
phenomenological, with many parameter values obtained by a fitting process. Little can be 
interpreted thereby about the intrinsic nature of the target, and the associated theory cannot 
be made predictive for estimations of the unmeasured (or unmeasurable) cross sections as 
needed in many other studies. That is also the case with many applications of .R-matrix 
theory of scattering [lj], though at least one 2], for n- 12 C scattering, sought to specify 
the background scattering from an optical potential with partial widths of the resonances 
specified from (p-s-d) basis shell-model wave functions of the nuclei involved, i.e. 12,13 C 

However, a more promising multichannel coupling theory of the scattering of neutrons 
from nuclei has been developed by the Padova group 0,0,IEIE0@|; a theory based upon the 
work of Rawitscher and Delic (3, [lOj • The approach is predicated upon making finite-rank 
separable representations of realistic interaction potentials and the properties of scattering 
matrices for such Schrodinger interactions. 

In brief, the approach starts with auxiliary sturmian functions (Weinberg states) forming 
a basis set in the interaction region. Initially one chooses a solvable potential problem at 
any suitable fixed negative energy from which first generation sturmians can be specified 
in closed analytic form. Second generation sturmians built upon the putative interaction 
potential matrices for a multichannel scattering problem of interest then can be found as 
linear combinations of that first generation set. The essential expansion coefficients result 
from a matrix diagonalization process. 

The scheme enables expansions (usually truncated to finite rank for numerical applica- 
tion) of the chosen interaction potentials in terms of those second generation sturmians, 
each in the form of a sum of separable interactions. The analytic properties of the scatter- 
ing matrix from a separable Schrodinger potential gives the means by which a full algebraic 
solution of the multichannel scattering problem can be realized. Of note is that the alge- 
braic structure of the Green functions, that lie at the heart of the S matrices, facilitates the 
identification of all resonances. The spin-parity (J n ), centroid, and width of each resonance 
can be ascertained without the need of a super-fine grid of energy values. It is important to 
identify especially narrow resonances that the scattering model predicts. 
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The multichannel formalism is outlined next in Sec. |Hj In this section there are two 
subsections, in the first of which we give some details of separable representations of potential 
matrices in terms of sturmian functions. In the second subsection, the process by which 
we can identify and locate resonances is described. The potential matrices specified using 
the Tamura collective model of scattering jll| then are discussed in Sec. 11111 This section 
also outlines how we deal with the Pauli principle, using orthogonalizing pseudo-potentials. 
Results of calculations in which deformation of that collective model for interaction potential 
matrices has been taken to second order in the case of neutron scattering from 12 C allowing 
coupling to the 0^ (ground), 2f (4.4389 MeV), and (7.64 MeV) states, and found for 
neutron (laboratory) energies to 5 MeV, are presented and discussed in Sec. IIVI Note 
that we use laboratory frame energies when dealing with the scattering domain as the data 
to be used is so specified. However, when considering the bound states (below the n- 12 C 
threshold) and/or spectra of 13 C we use energies in the center of mass system as such are 
used in tabulations of that spectra [l2T ]. 

But the collective model of nuclear structure is limited in scope and eventually we shall 
use this algebraic approach with microscopic (shell) model specifications of the target nuclei 
providing structure information which, when folded with an effective two-nucleon force, 
will determine the base input nucleon-nucleus potential matrices. That such shell model 
structures are reasonable for this purpose, at least with light mass targets, is considered 
finally in Sect. 



II. THE MULTICHANNEL T MATRIX FROM SEPARABLE INTERACTIONS 

Consider a system of C channels for each allowed scattering spin-parity J 71 " with the in- 
dex c (= 1, C) denoting the quantum numbers that identify each channel uniquely. Suppose 
c = 1 designates the elastic channel. The integral equation approach in momentum space for 
potential matrices V^Qo, g), requires solution of coupled Lippmann-Schwinger (LS) equa- 
tions giving a multichannel T matrix of the form 



°P en roo a;2 

£ L v cc>> (p> x ) p — zrr - T c>>c q;E)dx 

losed „oo ^ X 2 tt 

£ / Y cc» (P, x) 2 „ T£ d (x, q; E) dx 
,77^1 Jo hi, + x l 



c"=l 



Therein the open and closed channels contributions have been separated with the respective 
channel wave numbers being 



k c = \J ^{E - e c ) h c = ^/j(e c - E) , (2) 

for E > e c and E < e c respectively with e c being the threshold energy of channel c. Here \x 
designates 2m re d/fr 2 with m re d being the reduced mass. With the J n superscript understood 
from now on, solutions of Eq. (PJ are sought using expansions of the potential matrix elements 
in (finite) sums of energy-independent separable terms, 

N 

Vccf(p, q) ~ J2 Xcn(p) Vn 1 Xc'n(q) ■ (3) 
n=l 
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The link between the multichannel T matrix and the scattering matrix is 0, H| 



Sec' = bed - in^y k c k c ' T c 



N 

See' - 2 /c '^ c+1 7T/i V J k c Xcn (K) ( [V ~ G ] _1 ) , Xc'n' {K') J K> , 



(4) 



n,n'=l 



where now c, c' refer to open channels only. In this representation, Go and r\ have matrix 
elements (for each value of J n being understood) 



[G 



[rj 



"open 

E 

.c=l 



.r 



closed 



fc;? — x 2 + ie 



Xen' 



dx — / 



Xcn\Xj 



X 



hl + x 2 



Xcn'(x) dx 



(5) 



The bound states of the compound system are defined by the zeros of the matrix determinant 
when the energy is E < and so link to the zeros of {\rj — Gq\} when all channels in Eq. (J3J) 
are closed. 



A. Sturmian expansion of a multichannel interaction 



It is convenient to define first generation sturmians as solutions of uncoupled equations 
involving a known (local hermitian) potential matrix UjP\ 



(0) 



(0) 
-Vei 



(0) 



(6) 



where G^(E C ) is the free Green function evaluated at any suitable (arbitrary) negative 
energy E c . In these studies we chose 



E r = B 



-lMeV 



(7) 



which is independent of the channel. Then, as both operators G^ and UjP^ are hermitian, 



(o) 



U^Gf\E c 



-r/ 0) ($ {0) 

'lei \ ^ CI 



These eigenfunctions satisfy a potential orthonormality relation, 



(o) 



$(0)\ = Wg 

^C7 / 'ICl U - 



(0); 



(8) 



(9) 



where the normalization has been chosen to be consistent with a potential completeness 
relation of 



Ef/ (0) 



(0) 



(0) 



-1 



(0) 



Up) = up) 



(10) 



i=i 



A negative sign appears in Eqs. (jHJ) and (JHJ) so that the orthonormality and completeness 
relations are consistent with the convention that $2 ) are purely real functions. Details of 
how first generation sturmians may be evaluated have been published but are given in 
brief in Appendix El for completeness. 

Assuming that the infinite sums in the expansions Eq. (jlOj) can be truncated at a num- 
ber Ni of terms, sufficiently large that all important elements of the actual multichannel 
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scattering potential matrices can be expressed equivalently by either of two separable ap- 
proximations, one finds 



i=l 



(0) 



(0) 



(0) 



Ni 

V cc i ~ E Vcc' 
i=l 



c'l 



„(0) 
'Ic'i 



1 '*<°> 



= VK 



(i) 



(11) 



The procedure is to introduce second generation sturmians as eigenvectors of the coupled- 
channel homogeneous equations involving this first generation approximation , 



E gPv£> 



c'=l 



dp 



$(i) 



cp 



(12) 



where now it is assumed that the channel numbers are finite in extent (r). With the right side 
form of the first generation approximation for vj d , the expansion of the second generation 
basis in terms of the first is 

Ni 
3=1 



$ (0) \ 
cj / ' 



(13) 



where the coefficients are given by 



(14) 



These coefficients may be determined as solutions of a matrix equation that is formed by 
projecting the second generation sturmians given in Eq. (|T3*j) onto i^^ln V<j c and summing 
over the channel index c to find 



r iVi 

^ t ^ t ^c'm,cj Qcj,p T]p Qdm,p > 
c=l j=l 



where the cu-matrix elements are 



^dm,cj 



(^l^ C |<). 



(15) 



(16) 



If the channel coupling problem is assumed to be fully described by the selected set of V 
channels involved, the potential matrix V cc > is hermitian. The diagonalizing matrix Q then 
can be chosen orthogonal, with the T x Ni eigenvalues r/^ being purely real quantities. 
Then, a completeness relation 0, 0, E| for the potential V cc > can be established, namely 



rxTVi r r 

= E E EK C "|$S>KT & 

p=l c "=l c"'=l 

rxNi i 

= E \Xcp) — (Xc'p\ = Vcd , 
p=l T]p 



(1) 
d"p 



Vd"d 



(17) 



which provides the separable expansion of the potential matrix in terms of second generation 
(coupled-channel) sturmians. The form factors defined by 



hap) = E ^d $S> 



(18) 



c'=l 



6 



have been discussed in detail elsewhere |4|, |5|, 

With $ cra (r), the Weinberg (or sturmian) functions in coordinate space, those form factors 
then are given by 

Xc P (r) = £ f° V cd {rr')^ p {r>) dr' , (19) 

c'=l J0 

for a nonlocal potential matrix V cc r(r'r), while for a local interaction they have the form 

Xcn{r) = jzV ccl {r)^l{r) ■ (20) 

c'=l 

The form factors used to define the T and S matrices in Eqs. and (j3J) however are in 
momentum space. Those momentum space form factors are the Fourier-Bessel transforms 



Xcn{p) 



"2" 


2 * [ 


_7T_ 


p Jo 



Fe(pr)Xcn(r) dr , (21) 



where I is the orbital angular momentum quantum number. 

In application, it may be feasible to reduce the dimensionality of the problem by trun- 
cating the above to have N 2 < (T x Ni). The retained terms still are combinations of all 
of the r x Ni first generation sturmians however. Thus all relevant features of the fuller 
calculation can be included in a smaller basis calculation. The choice may reside with the 
listing of the eigenvalues in an ordered (decreasing magnitude) set. Also, note that the 
extension to higher generation expansions has been investigated 0, El El , but no significant 
improvement in the form developed with the second generation method was found. 

Formally, a set of coupled equations can be replaced by one for just the elastic channel 
alone from which one can define, in configuration space, the optical potential for elastic 
scattering. Of note is that, even assuming a local form for the elastic channel element of the 
potential matrix, the resulting optical potential will be energy dependent and nonlocal 0, 
El El EH- For completeness, a brief development of the optical potential resulting from the 
sturmian expansion method is given in Appendix [Bj 



B. Resonance identification 



In this section we describe a numerical technique for a rapid determination of all narrow 
resonances (in nucleon scattering from spin zero targets) arising from a system of coupled- 
channel Schrodinger equations. We consider only the elastic scattering channel for which the 
scattering matrix is recast (for each J n and with k — k±) using trivial matrix manipulation 
ofEq.(HJ). 



N 



S u = 1 - Z7r/i k Xin(k) (rj-G 



,Xln'(k) 



nn—1 
N 



mji k xin(k)- 



l-ry 2 Go?? 



nn'=l V 7 ?*"' 

Here, the diagonal (complex) matrix r\~ 2 is defined as 



n -l 
2 



1 



nn' \/f]n' 



--Xln'(k) 



(22) 



_ 1 

77 2 



(23) 
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The complex-symmetric matrix r/zGo^ 2 is then diagonalized 

N i ~ 

^ r] n "2 [G ] nn , Vn'~*Qn'r = CrQnr , (24) 
n'=l 

and the evolution of the complex eigenvalues ( r with respect to energy define resonance 
attributes. Resonant behavior occurs when one of the complex ( r eigenvalues passes close 
to the point (1,0) in the Gauss plane. From Eq. (}2*2*j) it is evident that the S-matrix has a 
pole structure at the corresponding energy where one of these eigenvalues approach unity, 
since one can write 

i-H N 1 

(l-rj^G r)^) = Qnr : _ T Qn'r • (25) 

J nn' r=l ' r 

These we designate as the resonance identifier equations in the text to follow. 

The eigenvalues ( r correspond to the positive-energy eigenvalues of the homogeneous 
(sturmian) problem with potential V^/(p, q) as given by Eqs. Q and ()17|). In operator form 
the relevant sturmian equations are now 

Y.Gf\E)v!$\$l{E)) = ( r (E)\€\E)) . (26) 

c' 

Note that here the sign convention is different with respect to Eqs. (JOJ) and (JHJ) so that the 
attractive eigenvalues will range in the upper-half Gauss plane as is conventionally used for 
resonance identifications. 

To show that ( r (E) and ( r (E) are the same, multiply both sides of Eq. (|2*S|l by the 
potential V^y , sum over the relevant channels, and use Eq. (JT7|) for V^J . The result is 

E IXcn)- [G (E)]nn< — (E(Xc'nWr(E)}) = ( r (E) £ \ Xcn ) - (^(Xcnl^ (E)}) . 
n,n'=l W W \ c ' / n=l Vn \ c ) 

(27) 

Projection onto the bi-orthogonal states then yields 

E ^ [Go(^)] n „, J- fE(Xc'"'l^(^))) = fE(XCn|0S(^))l , (28) 

n'=l »/« ^' V C ' / Vn\ c , ) 

which is equivalent to Eq. (|24|) provided that one make the identifications 

( r (E) = (r(E) ; Q m = ^= ( E<Xc'n| ) . (29) 

\/Vn \ d ) 

Thus the eigenvalues ( r of Eq. (}2*4*j) are the positive-energy sturmian eigenvalues of the 
potential V^)(p,q). Their general properties are well known (see, e.g., Ref. (Til) and they 
can be evaluated reliably (l5J. 



III. A COLLECTIVE MODEL OF THE POTENTIAL MATRICES 

The basic scattering model used was that for collective excitations as defined by 
Tamura llll . Such has been used before with some success to study resonance scatter- 
ing Il7| . and thus is our choice for this study. However, we allow certain extensions to 
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the usual collective model specification. Specifically, all terms to order '/3 2 ' in deformation 
have been carried given that the collectivity of the nucleus studied is strong. Also, the n- 
12 C potential field was allowed to have central (Vb), £ 2 -dependent (Vu), and spin-orbit (Vt s ) 
components. And, as suggested in Ref. [Hf, a spin-spin term (V ss ) has also been included 
in the central components. 

The basis of channel states is defined by the coupling 



\ C ) = \{t 8 )jir) 



\{£®s)j)® 4 a) 



M,7T 



J 



(30) 



where consideration is restricted to s = |, and the target states are denoted by ipi^ ) K a 
denotes any other quantum numbers necessary to uniquely define the target states). In 
coordinate space, the \{l®s)j) are the spin-angle functions of the relative motion wave 
function. 



A. The channel-coupling potential matrix 



Again with each J 71 " hereafter understood, and by disregarding deformation temporarily, 
the (nucleon-nucleus) potential matrices may be written 



VrjJr) 



{{£'s)J'I' | W(r) | (£s)JI) 

f{r)\v 5 c , c + V n [£ ■ £\ dc + V ss [s ■ I] c , c j + g(r)V ls {£ ■ s] 



(31) 



in which local form factors have been assumed. This potential matrix form results when one 
considers the basic Tamura collective model typically with Woods-Saxon form factors, 



f(r) 



n -1 



1 + e^ 



9{r) 



ldf(r) 



r dr 



(32) 



Deformation then is included by means of the rotational model 11] approach in which the 
nuclear surface is defined by 



R = Ra(l + e) 



E 

L(>2) 



47T 



2L 



Pl [Y L (f)-Y L (f )] 



(33) 



where T designates internal target coordinates. Expanding f(r) to order e 2 gives 



/(r) = /o(r)+e 



df(r) 
de 



d 2 f(r) 
de 2 



(34) 



There is a similar equation for g(r). Both f(r) and g(r) are now radial operators. 

In what follows that extended (collective) model is developed assuming that only one 
multipole L need be retained in Eq. (|33|); the definition of e. For the example we give later, 
of low energy neutron scattering from 12 C, that value is 2. However it is useful to specify 
the potential matrices for general L, as the results form a basis from which to develop the 
theory when more than one multipole is needed. 
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Using the property of tensor products [Hi ]. 

1 2L (or A- 1) 2 

[Y L (r)-Y L (f )] 2 = — E < L0L0 I £ °) 2 [^( f )- Y ^( T )] > C 35 ) 

£ even 

the radial operators are given by 



1 2L -I j2 x / \ 

+ -^(2L + 1) J; T^TT) (^OLOKO) 2 [Y r Y,] -Mi . (36) 

A similar equation applies to <?(r). 

To determine the deformed channel potential, it is not simply a matter of taking the ma- 
trix elements of the radial operators between channels states c and c' and substituting them 
into Eq. (pHj) . The channel potential expression involves matrix elements of the products of 
two operators and so one must first make symmetric the potential matrix form and use a 
suitable completeness relations of the type J2 C \ c >< c| = 1 to separate the action of product 
operators. The correct form of the matrix V c / C (r) is 

V c , c (r) = V f ctc (r) + ^ ii E{[«]c'c»/c» c (r) + / cV (r)[«] c " c } 

1 c" 

+ o Wis E {[^ s ]cc"9c"c(r) + gc'c"(r)[£-s] c/ > c } 

Z c" 

+ ^EiMW/c-cM + /cV'MMU} , (37) 

1 c" 

where W\ a = 2VJ s A 2 r . Using the Woods-Saxon forms, details of these matrix elements are 
given in Appendix O 

A feature of our calculations is that the deformation is taken in all terms of the potential. 
In Ref. (llj . only the central part of the potential was taken to be deformed. This is an 
important point. It is incorrect to disregard the effect of deformation on the non-central 
parts of the potential, particularly so if large deformation values are to be used. 



B. The Pauli exclusion principle 

The standard coupled-channel approach to nucleon-nucleus scattering involves elastic and 
inelastic processes between two particles where at least one of the two is itself composite 
and interacts with the other through effective (optical) potentials. One must keep in mind, 
however, that the underlying process is a complicated many-body scattering problem which 
requires consideration of the Pauli exclusion principle. 

It is a known fact that phenomenological (coupled-channel) models coupling collective 
deformations with single-particle optical potentials violate the Pauli exclusion principle j2(j . 
In single-channel scattering processes this violation does not represent a severe problem 
because the exclusion principle can be taken into account in the scattering process naturally 
if the single-particle potential contains a series of deep bound states (the forbidden states) 
to which the scattering wave function is orthogonal by construction. 
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In a coupled-channel model another problem arises because, in addition to the orthogo- 
nality of the scattering wave function to a series of forbidden states in the elastic channel, 
one has to eliminate also all virtual transitions to the forbidden states in the excited chan- 
nels. A failure to do so implies the construction of an erroneous scattering wave function; 
one which would be unavoidably contaminated by the unphysical couplings to states of the 
nucleon- (excited nucleus) system forbidden by the exclusion principle. 

The main effect of the Pauli principle is well represented by the suppression of a part of 
the phase space which would be otherwise accessible to the system. This implies that the 
two fragments being scattered must negotiate between themselves with effective potentials 
where certain states forbidden by the Pauli principle are excluded explicitly. 

An efficient method to achieve this goal is obtained by introducing the orthogonalizing 
pseudo-potentials (OPP) [2l|,|22| The OPP method is a variant of the orthogonality condition 
model by Saito |23J] which allows projection of the scattering solution onto the subspace 
permitted by the exclusion principle through a suitable renormalization of the two-particle 
interaction matrices. All virtual transitions to the forbidden states are therefore eliminated 
by the redefinition of the two-body interaction operators. 

Implementing this method we add an orthogonalizing pseudo-potential (OPP) term to 
the collective potential matrices V cc '{t) giving 

V«,(r, r') = V cd (r) S(r - r') + XA c (r)A c (r f ) S c , c , . (38) 

The pseudo-potential is manifestly nonlocal in coordinate space and it is diagonal in the 
channel index. The function A c (r) is the radial part of the single particle wave function 
in channel c, spanning the phase-space excluded by the Pauli principle in that particular 
channel. In the actual calculation, A c (r) is determined by solving the radial Schrodinger 
equation for each channel c 

d 2 

— A c {r) + w c {r)A c {r) = , (39) 
numerically with bound-state boundary conditions. Here the abbreviation 

w c (r) = 2 J[E- V cc (r)} - . (40) 

has been used. The OPP is added only for those channels c containing single particle quan- 
tum numbers referring to closed shell configurations. For the n- 12 C system to be considered 
here, we eliminate all elastic and inelastic states related to the Osi and On 3 configurations 

2 2 

(deep bound states) from the scattering equations. Following the OPP method, the forbid- 
den configurations are eliminated from the dynamical equation in the limit A — > oo. For this 
study on low-energy nucleon-nucleus scattering, we have selected A ~ lOOMeV as a value 
sufficiently large. 

It is an important feature of the sturmian expansion method that one obtains an algebraic 
solution for coupled-channel scattering for both local and nonlocal interactions. The physical 
nonlocalities introduced by the Pauli principle are reflected in the strong nonlocal character 
of the OPP terms. Thus, we apply the sturmian expansion method to the Pauli corrected 
potential V cc '(r, r') and not to the original potential V cc /(r). This leads to the following 
modifications of the sturmian expansion method: 

1) The matrix uo c i m c j now has to be calculated according to 

fVmrf = (^l |Vcc| , (41) 
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however, the auxiliary sturmian base Q c i m j can be maintained as defined before. 

2) Diagonalizing the matrix u leads to the coupled sturmian eigenstates ^p ) an d eigen- 
values rffl, thus providing the ingredients for the expansions of the potential matrices 



rxTVi r r 

= E E E v cc » \*%) 

p=l c"=l c"'=l 

E \Xcp) — (jy (Xc'pl j 



p=i f?p 
where the last equation implies 



(i) 

e"'p 



v, 



v. 



(2) 



(42) 



c=l 



(43) 



Therefore, the a;- matrix is now calculated with the integrations 



^c'm,cj 



^i(r)V c , c (r)^>(r)dr 



(0), 



A c (r)$W(r)rfr 



A c (r)$Jf(r)dr 



(44) 



and the new potential form factors in coordinate space are given by 



X»(r) = £Mr)*£>(r) + AA 



c'=l 



A c (r>^V)dr' 



(45) 



IV. APPLICATION: n- 12 C SCATTERING 



As a test case for study using the algebraic approach with the collective model potential 
matrices, we consider neutron scattering from 12 C. The measured cross section data are 
shown in Fig. ^ Those values were obtained using CINDA search in the web page of the 
National Nuclear Data Center, Brookhaven (www.nndc.bnl.gov). References for all the data 
sources are given therein. The data set shown reveals very narrow and very broad resonances 
in the energy range from to 4 MeV. Most prominent of the narrow resonances are those at 
2.08 and at 2.75 MeV. They have been assigned spin-parities of | + and | + respectively. The 

broad peaks have spin-parity assignments of | + , | + , and \ centered at 2.95, 3.58, and 4.26 
MeV respectively. Those resonances lie on a background that varies smoothly with energy 
from 4.6 barn. 

In our analyses, three states in 12 C have been taken as active. They are the 0^ (ground), 
the 2f (4.4389 MeV), and the 0^ (7.96 MeV) states. As the three states have the same (+) 
parity, the even and odd parity scattering channels then depend separately upon the even and 
odd parity input potentials respectively. To begin the search for suitable collective-model 
parameters, we made the usual assumptions of the shell model, and nuclear deformations; 
namely, that the nuclear radius is given by Ro = roA 1 ^ 3 , taking ro = 1.35 fm, the diffuseness 
as a = 0.65 fm, and that the depth of the potential be about Vq = —40.0 MeV. The 
strength of the spin-orbit potential was taken as V^ +) = 6.7 MeV. The value (3 2 = -0.6 
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FIG. 1: Measured elastic cross sections for n- 12 C scattering as functions of neutron energy. 



was taken to define the deformation of the 12 C target nucleus. As usual, the incoming 
nucleon is treated as a point particle. The other parameters, the strengths of the £ 2 and 
spin-spin interactions, are not well known, and were used as fitting parameters. As well, the 
other potential strengths and the deformation parameter were considered as parameters that 
could be varied to improve the representation of the data. A first crude attempt to represent 
the experimental results led us to increase the positive-parity central-potential strength to 
V + = —46.848 MeV, set the £ 2 term strength at = 0.611 MeV and both spin-spin term 
strengths at = — 1.0 MeV. The first generation sturmians were obtained using a square 
well with parameters {B, R} = { — 1.0 MeV, 7.0 fm} as the auxiliary potential. Matrix sizes 
were limited to a 30 sturmian expansion for each channel. The resonance identifier equations 
as well as those for the S matrices were evaluated for all channel spin-parities J n from | 
to | . Using this starting set of parameter values, by taking deformation through second 
order, lead to a very rich structure in the scattering and more importantly a structure that 
is reflected in measured data. 

The scattering model we use does indeed have the features seen in measured data though 
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the centroids and widths can, and do, vary considerably with choice of parameter values. 
Also, dependent upon the starting parameter set, there can be more narrow and/or weak 
resonances in this energy regime. However it is essential that a complete coupled- channel 
analysis be made of the scattering for the resonances to be considered sensible. Of course 
the deformation parameter (3 2 has been assessed from diverse applications of the collective 
model of excitation |24[ and in particular from the electron form factor and so should be 
rather constrained in any variation. Likewise proton inelastic scattering cross sections from 
excitation of the 2 + (4.4389 MeV) state indicate potential parameter values that one may 
consider as "sensible". But those analyses, by and large, used the (distorted wave) Born 
approximation. Such an approximation might equally well be used in finding the background 
at low energies but is not appropriate when a study of the resonance attributes are to be 
made. 

The end result of the search process we have made is the set of parameter values 

= -49.1437 MeV 
Vlf ] = 4.5588 MeV 

= 7.3836 MeV 
VH = -4.77 MeV 
ro = 1.35 fm 
a = 0.65 fm 

Initially the search process was extremely computer time consuming since, prior to the 
development (and use) of the resonance identifier equations, either the phases of the S(± 
(elastic) scattering matrices and/or the calculated cross section for many energies to 5 MeV 
were used to specify the resonance energies and their FWHM. But using the calculated 
cross section and/or the phase properties of the elastic channel S matrices alone does not 
guarantee that weak and/or very narrow resonances will be evident. The choice of energies 
and the step size used may not reveal characteristic effects that draw one's attention to the 
energy region for more intense study. It is a hit or miss scenario. Indeed with the initial 
coarse grid of energy steps of 0.01 MeV even the dominant | + resonance near 2.1 MeV MeV 
was not evident in the base calculations that were made. If there were no other means by 
which resonance existence and centroid energies generally could be located, predictions for 
poorly or as yet to be measured cross sections would have to be made using inordinately 
small energy steps over the whole range of interest; a major computing problem. Fortunately, 
the properties of the Green functions and of their eigenvalues enable the existence, number, 
and energy centroids of resonances to be found at the outset. 

Using the resonance identifier equations with sturmians built from the basic interaction 
potential, we obtained a sequence of bound states and resonance energies and widths for 
resonances in the neutron plus 12 C system. The actual values and how they are specified 
are presented in the second of the following subsections. In the first we illustrate the general 
behavior of the sturmian eigenvalues in Argand diagrams where the horizontal axis gives 
the real part of the eigenvalue and the vertical, the corresponding imaginary part. 

1. Sturmian trajectories 

For each total angular momentum and parity, we have calculated the eigenvalues ( r {E) 
over the energy range 0.01 to 4.96 MeV with a constant step of 0.05 MeV. The eigenvalues 
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-47.5627 MeV 
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V£~> = 9.1760 MeV 
y(+) = -0.052 MeV 



(46) 



3 2 = -0.52 . 
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FIG. 2: Argand diagrams for the energy variation of the sturmian eigenvalues for J 71 " = l/2 + . 

in this energy regime, usually increase in magnitude with energy, with those coinciding with 
definable resonances in cross sections having values less than 1 + iO for energies below each 
associated resonance energy. Thus even starting with a fairly coarse energy grid, as long as 
there are energy values both below and above the resonance centroid, no matter how narrow 
the resonance, an eigenvalue will change from below to above 1 + i0 for the two energy points 
in the grid that lie below and above that centroid. Thus we learn from use of the resonance 
identifier equations not only how many resonances there are in the energy range considered, 
but also where to make finer grid searches to better ascertain the characteristic properties 
of each resonance. The results to be shown were all derived from the calculations made with 
the potential matrices associated with the parameter set of Eq. 0(3 

In Fig. |21 we show the behavior of the four largest eigenvalues for J n = | + in the energy 
range to 5 MeV. The unit circle is displayed by the dashed curve. Note that the graph is semi- 
logarithmic since all of these eigenvalues have small imaginary components. That is due to 
the elastic channel always being open. Some trajectories continue below the graphed range. 
All four at low energies have a variation that is vertical to the real axis; a characteristic 
of the eigenfunction solutions for s-waves. The trajectories of two eigenvalues (curves (1) 
and (4)) show typical behavior of single- channel (potential) s-wave eigenvalues, as has been 



discussed in Ref. [14j. They do not correspond to any resonant structure. The trajectories 
of the other two eigenvalues (curves (2) and (3)) shown in Fig. El exhibit a quite different 
behavior; one which originates from the coupled channel dynamics. In the energy range 
considered only that identified as (2) links to a resonance feature in scattering. 

Next, the behavior of J 71 " = | eigenvalues are considered. The relevant Argand diagrams 
are given in Fig. El The point values of the largest eigenvalue (labelled (1)) has been i 
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FIG. 3: Argand diagram for the energy variation of the sturmian eigenvalues for J 71 " = 3/2 + . 

connected by a long-dashed line to guide the eye to link the energy sequence of the results. 
The actual trajectory is a cusp. That is very evident with curve (1), but similar though 
more slight features are evident in the trajectories (2) and (3). That cusp feature of the 
eigenvalue trajectories links to the opening of the 2f state at 4.4389 MeV in the coupled- 
channel algebra. 

Note that these eigenvalues again have small imaginary parts and so the plot is semi- 
logarithmic and again some details of the trajectories lie below the graphed range. In this 
case the trajectories do not depart vertically from the real axis at the scattering threshold 
since they are <i-wave solutions. The (largest) eigenvalue clearly evolves well beyond the unit 
after crossing with a quite small imaginary part. Thus it coincides with the lower energy | 
compound resonance in the cross section. Likewise the 2 nd trajectory crosses the unit circle 
at higher energy and with a larger imaginary part. This coincides with the known second, 
broader, | + shape resonance at 3.4 MeV. The 3 rd and 4 th sturmian trajectories shown in 
Fig. El track towards the unit circle but have not crossed before 5 MeV. 

Finally, in Fig. EJ we show the Argand diagram of the 2 nd and 3 rd largest eigenvalues 
for J 71 " = | + since the first coincides with a state in the mass 13 spectrum below threshold. 
The unit circle is represented therein by the dashed curve. Again the long-dashed line is 
simply to guide the eye to the energy variation of these results; results which give rise to 
a very narrow resonance at E~2.086 MeV. As the vertical scale again is logarithmic, the 
trajectories continue below the graphed limit. The trajectories also exhibit a non-vertical 
"take-off' from the real axis commensurate with them portraying d-wave eigenvalues. Also 
both exhibit cusps which tag to the energy of the 2f state. 

Given that the imaginary component of the eigenvalue (curve (2)) is very small, it is a very 
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FIG. 4: Argand diagram for the energy variation of the sturmian eigenvalues for J 77 = 5/2 + . 

narrow resonance in the cross section, and the resonance centroid essentially is equivalent to 
the energy at which the real component of the eigenvalue itself is unity. That resonance in 
the total cross section is shown in Fig. |SJ It has a width of about 15 keV and a magnitude 
of over 6 barn. 



2. Calculation of resonance parameters 

To find the resonance centroids and widths, one first has to find all the energies where 
the real part of the sturmian eigenvalues crosses unity. That determines an approximate 
value for the resonance energy Er by the condition 

( r (E R )=l + i8. (47) 

Since these eigenvalues change smoothly with energy, it is possible to find these points with 
a rapidly converging predict-and-correct iterative procedure. We have also estimated the 
width of these resonances, according to the approximation 

- 2 x wm 1 (48) 

where d[R(( r )] represents the differential of the real part of ( r at the resonance energy. 
These formulas (Eqs. (}4Tj) and (|48j)) are correct in the limit 5 << 1, therefore they should 
be considered reliable only for the narrow resonances. For the case of wider resonances, it is 
not difficult to find the more complicated expressions required. In the context of the present 
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FIG. 5: The J 77 = | + resonance in the elastic scattering cross section. 

discussion, however, they are not needed and indeed cross-section evaluations at reasonable 
energy steps can be used to ascertain them quite easily. 



A. The theoretical elastic cross section and polarizations 

The method described in the previous subsection enables us to determine with some pre- 
cision, the centroid energies and widths of all the resonances. Similarly, we can find the 
energies of the bound states. At the same time we have sought a good theoretical descrip- 
tion of the elastic cross section in compare with the data of Fig. ^ The polarizations (at 
select scattering angles) will follow without having played any part in determining potential 
parameter values and so are "predictions" . 

Starting from the basic potential parameters, we determined all the resonance energies, 
widths and bound-state energies. That set we denote as Then we changed each 

parameter by a small amount and calculated a matrix of (approximate) partial derivatives 
of the with respect to the parameters. That matrix of derivatives we denote generically 

d p >. Then, to first order, a new set of energies |.E^ j are given by the n — 1 form 
of the recursive formula, 

ap(n-l) 

E (n) „ + ^ A(3p (4g) 

p (ypp 

where Aj3„ are the changes in parameters f3 p needed to produce the new E\ n \ If the number 
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of parameters is equal to the number of energy values, this set of equations can be considered 
as a linear system to solve for the A[3 P needed to produce the resultant E^ 3,1 '. We chose 
a few of the parameters to vary, and a few essential features of the resonance structure, as 
given in the data or read off from Fig. Q in an attempt to produce a fit to that data. 

This has turned out to be a rather complex task because of the interplay among the 
parameters, as well as due to some lack of flexibility of the collective model on which this 
analysis is based. It can be compared to squeezing a balloon . . . something always pops out 
that one wants to keep in. In fact, we carried out this process in two stages. In the first, a 
quite rough agreement with the data in Fig. [T]was found. This set of parameters was used as 
a "new base" and the process was repeated. In the end, we obtain a quite good description 
of the experimental cross section. But some small defects remain, so this result may not be 
the best that can be achieved. With this mult i- dimensional non-linear problem, there may 
be a number of "quite good" sets of parameters; we have found one. 

Rather than attempting to fit everything, we decided to limit ourselves to trying to 
reproduce the most prominent features of the experimental elastic cross section shown in 
Fig. [H and of the experimentally best established resonances. Specifically, we focused on the 
two | + resonances in the range 2.5 to 4.0 MeV energy, the prominent narrow | + resonance 
just above 2.0 MeV, and the known bound states below threshold. The collective model 
used produces two | + resonances in agreement with the data. Of the two, one is very wide 
(r ~ 1 MeV), likely single-particle, and can account for the broad peak centered around 3.5 
MeV. The other is narrow, generally less than 100 keV, and by interference with the first, 
produces the prominent structure in the cross section near 3.0 MeV. This is shown by the 
continuous curve in Fig. El In this region a narrow 1 resonance also is obtained from the 
model. Such is not seen in the data nor does a partner state exist in the spectrum of 13 C. 
Note that in Fig. El we use the evaluated cross section data file (ENDF) also found in the 
website of the National Nuclear Data Center, Brookhaven (www.nndc.bnl.gov). 

Nevertheless, and as can be seen in Fig. El the overall agreement between experimental 
cross-section data and theory is good. The theoretical curve shows the structure just below 
3.0 MeV as very similar in shape. The prominent narrow | + is reproduced at the correct 

energy. The smaller | + is observed in this figure too, and it is near the correct position. 
At low energy above threshold, the theoretical cross section is in good agreement with the 
smooth background curve seen in the experimental data. The theoretical curve also displays 
the bump above 4.0 MeV, though it is higher and broader than the data. This resonance is 
assessed to have spin-parity of | . A more careful study of the negative-parity resonances 
at higher energy is required. An appropriate treatment of higher-energy negative-parity 
states in the collective model should include the 3~ state at 9.641 MeV in our selected set of 
channels. By so doing, the dimensionality of the problem increases and computing time with 
the existing code becomes exorbitant. A new code based on massive parallel architecture is 
required. Currently such is under construction. 

Polarization data (at two scattering angles) are shown in Fig. They are compared with 
results that come from the S matrices determined by just seeking cross section features. 
The known resonances and a characteristic change of sign through the broad | + resonance 
centered near 3.5 MeV are well replicated, as is the variation attributable to that centered 
near 3.0 MeV. We stress that these polarization data were not used in seeking a best set 
of potential parameters and were the result of but a single calculation (for each scattering 
angle) using the cross section defined S'-matrices. 

The full set of resonance energies and widths and bound-state energies between -35.0 MeV 
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FIG. 6: Comparisons between theory (dotted curve) and experimental data (solid curve, ENDF 
evaluated) of the cross section from the elastic scattering of neutrons from 12 C. 
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FIG. 7: Comparisons between theory (solid curves) and experimental data (filled circles and 
squares with error bars) of the polarizations from the elastic scattering of neutrons from 12 C. 
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TABLE I: Resonances and bound states found using the resonance identifier equations and the 
parameter set of Eq. (|46|) . 
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and +5.0 MeV obtained with these parameters are given in Table |l] It is worth stressing 
that only when Pauli blocking is treated by including the OPP does the approach give four 
states of the correct spin-parity and at appropriate energies below the n+ 12 C threshold. 
Results shown in Table H] for the spectra in the continuum are those obtained using the OPP 
terms and are also those displayed in Fig. El 

The last two columns in Table gives the results obtained when Pauli blocking is ignored 
(by omitting the effects of the OPP). Those results are identified also by the header, "without 
OPP". The bound state spectrum that results most clearly shows the spuriosity. Far more 
states result than are known empirically (for 13 C), and most are more deeply bound than 
the known ground state. With the scattering cross sections, that spuriosity is not as well 
displayed. Indeed by adjusting parameter values we could attain nearly as good a result 
as found when the OPP is used. However, without the OPP the resonance states are not 
the correct entries in the coupled-channel sequence. For example, without the OPP, the 
narrow | resonance near 2.0 MeV is the third in sequence rather than the second when 
Pauli blocking is used. But, the high spin-parity resonances should not be, and are not, 
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affected by Pauli blocking. For example, in both calculated spectra the | resonance has the 
same centroid and width values. 

Focusing on the results with the OPP corrections, clearly the essential features of the even- 
parity resonances are reproduced well. However, this analysis suggests that the experimental 
resonance at 4.259 MeV would have a J n assignment of 5/2 + and that a broad resonance 
that does not influence the elastic scattering has a centroid at 4.833 MeV and would have a 
spin-parity of l/2 + . Experiment suggests that the l/2 + is observed and lies at 4.259 MeV 
while the 5/2 + is not to be seen in this energy regime. But it must be remembered that 
we have limited the process to just three states in 12 C, and have some allowance in the 
parameter specifications in the collective model prescription for the interaction potentials. 
Likewise there are two negative-parity resonant states indicated from these calculations; 
namely the l/2~ "hidden" in the structure above 3.0 MeV and a 5/2~ which might align 
with what is a suggestion of a small narrow resonance just above threshold in the data. 

Despite the few anomalies, it is clear that this method makes possible exploration of 
resonance behavior, in this first presentation by using the collective model, and by using 
purely algebraic means. However, it is also clear that the collective model, coupled with the 
restrictions imposed by limiting the choice of active target nucleus levels, will never provide 
as complete agreement with experimental data as one would like. It is planned to define 
interaction potential matrices from folding a suitable low energy NN force with microscopic 
model structure such as given by the shell model; a process that has had much success in 
recent years in analyses of higher energy nucleon-nucleus scattering |25|. Of course Pauli 
effects would be approximately included via this folding when full antisymmetrization is 
considered. Therefore, it is of interest at this stage to speculate what might be learned from 
a more microscopic theory, and that is discussed in the following section. 



V. THE STRUCTURE OF 13 C IN TERMS OF 12 C PLUS A NEUTRON 

The algebraic scattering program has been applied to low energy neutron scattering from 
12 C. The energy variation of the measured cross section reveals sharp resonances, some of 
which one may associate with the presence of states in 13 C. Likewise, the algebraic methods 
allow prediction of bound states of the n- 12 C system which should correspond to states in 
13 C below threshold. In our first calculation, the structure of 0^ (ground), 2f (4.4389 MeV) 
and 0^" (7.65 MeV) states of 12 C have been taken into account. That such a model should 
give results characterizing the observed cross section follows from consideration of (p-s-d) 
shell models for 12 ' 13 C. 

The key quantity that relates the structure of the spectrum of 13 C to the spectrum of 
12 C plus a neutron in a single particle orbit j is the pickup spectroscopic amplitude, 

S{ IJ = (¥j X) (1, 2, ■ ■ • , 13) ||a] J $,(1, 2, ■ • • , 12)) . (50) 

The development applies equally to a proton with the compound system then being 13 N. 
The particle coordinates will be omitted hereafter with the identification: 

\¥j X) ) = the A th state of 13 C with spin - parity J w 

|$/) = the state in 12 C with spin / . (51) 

Thus, while just three unique spin states in 12 C are considered, there will be more than 1 
state of spin J n formed by attaching any specific orbit nucleon. Whether these numbers 
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are predicted by a shell model calculation or instead are extracted from scattering data 
analysis, there are two sum rules that indicate if the coupling of any single particle state in 
the spectrum has been exhausted. 

The first is the pickup sum rule. It is the result of summing the spectroscopic probabilities 
over all target ( 12 C) states and is 
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On closure over the target states, this gives 
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where n*p is the number of nucleons of the appropriate orbit in the particular 13 C state 
\^ jm)- A further summation over all possible values of j gives 
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where N is the number of nucleons of the selected type (neutrons in the case considered). 
For the 13 C-[ 12 C + n] systems, spectroscopic amplitudes have been calculated jH using a 
shell model in which the active shells were Op-ls-Od. Recently, the structure of 12 C has 
been found [26] using a complete no-core (0 + 2)hu shell- model. Such has not been used 
to specify 13 C as yet. However, the basic details should not be too different to what is 
discussed. In S^-j 7 expressing the states of 13 C in terms of 12 C plus a neutron are given. 
They are presented in the order of excitation in 13 C with the excitation energy shown in 
column 2. The equivalent energy in the center of mass for n + 12 C is shown in column 3. 
The first four entries therefore are closed to neutron scattering from 12 C. In the incident 
energy regime to the threshold of excitation of the 2f, (4.4389 MeV) state in 12 C, four 
known resonances are expected coinciding with the spin-parities of the compound system 
being f + |(2); f + |(i); | + l(2) and \ |(2). The third of those however is very broad, O(MeV), 
and .R-matrix studies |2| suggest that it is a single particle potential resonance. 

Thus the shell model indicates that the four known states in 13 C below the n+ 12 C thresh- 
old are well represented by a neutron coupled to the ground and 2f states of 12 C and so 
should be (and are) well defined by the negative energy solutions of the multichannel alge- 
braic scattering problem. Likewise, coupling a neutron to the same states in 12 C give shell 
model candidates for states in 13 C above the n+ 12 C threshold; ones that match observations 
from scattering data. They relate to resonances in the scattering of neutrons from 12 C which 
the multichannel scattering theory also defines well. The sizeable components of a neutron 
coupled to the 2 + state of 12 C that the shell model indicates is consistent with the strong 
coupling we have found necessary via a collective model prescription in these evaluations of 
scattering observables. 
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The second sum rule, the stripping sum rule, is formed by completing a sum over all 
compound mass ( 13 C) states. It is important as it indicates whether or not the model 
calculation may give the large components of any state. That sum rule is defined by 
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(55) 



Using the (anti) commutation property of the creation/annihilation operators this sum re- 
duces to 

S 5 = (2J + l) 2 [2j + l-n i ] , (56) 

where rij is the number of nucleons in the orbital j in the target ( 12 C) state \<&ik)- in 
practical calculations one cannot deal with all of the mass 13 levels and so this sum rule is 
best viewed as an estimate of how much transition strength lies with states other than those 
investigated. 

Using the p-s-d shell model gives stripping sum rule values listed in Table UTTl They are 
compared against a theoretical limit set assuming that there are no Is — Og? neutrons in 
the ground state description of 12 C. That is quite reasonable as a shell model study [26| 
made using the complete (0 + 2)hu space gives 11.6 nucleons within the 0s — Op shells for 
both the + (gs) and 2 + (4.4389) states. However, with the simplest of shell models (packed 
orbits), the pickup sum rule values theoretically are 2, 4, and 6 for the lsi/ 2 , 0d 3 / 2 , and 0d 5 / 2 
shells respectively. The occupancies of the 0si/ 2 , 0p 3 / 2 , 0pi/ 2 used to get the results listed in 
Table fTTTI are those given by p-s-d shell model calculations 0, namely 2.0, 3.288 and 0.712 
respectively. For those single particle specifics the theoretical sum rule limits are 0.0, 0.712, 
and 1.288. In the 2 + (4.4389 MeV) state in 12 C, the occupancies are varied slightly from 
those with the Op shell numbers most altered. The relevant theoretical sum rules then were 
found using occupancies of 2.0, 3.01, and 0.98. The sum rule exhaustion suggests that the 
assumption that the states in the listings above are a nucleon plus the 12 C nucleus in either 
the ground or the 2 + (4.4389 MeV) state is quite well satisfied except if that nucleon is in 
the g?3 orbit. It would not surprise, therefore, if the excitation of the 13 C | + states may not 
be as well described as others with any model involving the ground and 2 + states of 12 C in 



the basis. Note also from Table HT] that the ^ 



7 + 



and 



| states anticipated to lie at 2.54 and 
2.6 MeV excitation in the n- 12 C system, are solely based upon couplings with the 12 C 2 + 
state. As such they may be only weakly excited in the n- 12 C scattering cross section. 



VI. CONCLUSIONS 

A model for nucleon-nucleus scattering that uses sturmian expansions of multichannel 
interactions between the colliding nuclei has been used and found to give resonance scattering 
upon an average background; typical of what is measured with low-energy experiments. That 
sturmian expansion theory also gives resonance identifier equations whose solutions identify 
feasibly, all possible resonance aspects to the scattering problem. Indeed we have found 
that it is quite practical to find the narrow resonances (compound and quasi-compound) 
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TABLE II: Shell model values for pickup spectroscopic amplitudes giving C for neutrons on C 
with energies to 4.0 MeV. Single particle states belong to the Op — Od — Is space. 



13 C 


Eea; 


E rel. CM. 


neutron 


12 C 


One channel 


Two channels 


J (X) 


in 13 C 


n + 12 C 




I 


qjOj 


b W 


1-| 


0.00 


-4.95 


Pl/2 
P3/2 



2 


1.127 


1.107 
-1.498 


1 + 1 

2 ll 


3.09 


-1.86 


•51/2 
^3/2 
^5/2 



2 
2 


1.349 


1.349 
-0.148 
0.069 


3~l 
2 ll 


3.68 


-1.26 


P3/2 
Pl/2 
P3/2 



2 
2 


-0.840 


-0.887 
-1.819 
1.122 


5 + | 
2 ll 


3.85 


-1.09 


d^/2 
Sl/2 
d 3/2 
d$/2 



2 
2 
2 


2.271 


2.271 
0.334 
-0.293 
0.415 


5 + | 

2 b 


6.86 


1.92 


d<r,/2 
Sl/2 
d 3/2 



2 
2 


0.409 


0.409 
-2.173 
0.234 


7 + 1 
2 ll 






d-5/2 


2 




0.502 


7.49 


2.54 


d 3/2 
d<r,/2 


2 
2 




0.778 
-2.565 


2 ll 


7.55 


2.60 


Pl/2 
P3/2 


2 
2 




1.314 
0.495 


3 + 1 
2 ll 


7.68 


2.73 


d-3/2 
Sl/2 
d 3/2 



2 
2 


0.944 


0.944 
1.638 
0.192 


3 + 1 
2 12 






db/2 


2 




0.326 


8.20 


3.30 


d 3/2 
Sl/2 
d 3/2 
d b/2 



2 
2 
2 


-1.000 


-1.000 
0.803 
0.318 

-0.644 


1-| 
2 12 


8.86 


3.91 


Pl/2 
P3/2 



2 


0.082 


-0.049 
0.084 


3~l 
2 12 


9.50 


4.55 


P3/2 
Pl/2 
P3/2 



2 
2 


0.052 


-0.075 
-0.102 
1.123 


3~l 
2 13 


9.90 


4.95 


P3/2 
Pl/2 
P3/2 



2 
2 




0.031 
0.058 
-0.205 



25 



TABLE III: The sum rules for the two channel case - neutron space Op — Is — Od. 



Neutron 


12 C 


Theoretical 


Shell model 


exhaustion 


state 


state 


sum rule 


sum rule 


% 


s l/2 


0+ 


2 


1.820 


91 




2+ 


10 


8.387 


84 


Pl/2 


0+ 


1.240 


1.239 


100 




2+ 


5.080 


5.050 


99 


P3/2 


0+ 


0.763 


0.758 


99 




2+ 


4.935 


4.871 


99 


d-3/2 


0+ 


4 


3.110 


78 




2+ 


20 


0.95 


5 


d-5/2 


0+ 


6 


5.557 


93 




2+ 


30 


25.281 


84 



of the coupled-channel Schrddinger problem by diagonalizing the energy-dependent matrix 
rj i G (E)rj~ 2 and studying the trajectories of the relevant eigenvalues in the Gauss plane. 
When one of these complex eigenvalues evolves past the point 1 + iO and does so having 
also a small imaginary component, one of these resonances occurs. Then Eqs. (j4"Tj) and 
(refg-res) can be used to determine the resonance parameters. Since these eigenvalues have 
smooth energy dependencies, it is generally much simpler so to determine the occurrence of 
resonance states than by other means that have been used to date. 

We have used a collective model (to second order) to define a multichannel potential 
matrix for low-energy neutron- 12 C scattering allowing coupling between the 0+ (ground), 
1\ (4.4389 MeV), and 0^ (7.64 MeV) states with coupling taken to second order. The 
algebraic S matrix for this system has been evaluated. Good results have been found for 
the sub-threshold bound states and for the cross sections and polarizations as functions of 
energy. The latter are rich in structure having both narrow and broad resonant features for 
different J n ; the existence and parameter values of which can be ascertained with ease. 

The introduction of the Pauli exclusion principle in the coupled-channel model by means 
of the orthogonalizing pseudo-potential (OPP) technique ensures that the model gives a 
spectrum (bound states and resonances) that are built as physical states. This was an 
important improvement giving an unique identification of states with respect to experimental 
data. 

The results indicate that this approach has predictive power and can be used to interpret 
the experimental resonance spectra of the nuclear processes at low energy. The discussion 
on the n- 12 C resonant spectrum contained herein is just one simple example of the potential 
use of this approach in nuclear analyses. 
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APPENDIX A: THE FIRST GENERATION STURMIANS 



With I being the orbital angular momentum quantum number that is embedded in the 
set collectively denoted by c, the eigenvalues qu of the first generation sturmian problem 
may be calculated as follows. Consider a square well problem for a binding energy B and 
well with depth Vo and radius R. With xu = qup, q = \fpB, and y = qp, the first eigenvalues 
are defined from solution of the transcendental equations, 



xu cot xu + y = 



and 



xu- 



Fi-i(xu 



+ y 



; I > . 



(Al) 



(A2) 



F^xu) * n t (y) 

The functions Fi and Qi are the Riccati-Bessel and modified Riccati-Bessel (of the 3 rd kind) 
functions respectively; each of order I. They can be solved using well known recursion 
formulas 3]. In terms of these, the first generation sturmians are 



The normalization constant is 



Fi(qi n r) 



Fi (gup) 



Qi{qr) 



for < r < p 
for r > p . 



(A3) 



A, 



and the sturmian eigenvalue is 



{ql + q 2 )tfF?{q u r)dr ' 
pV 



(0) 
Vci 



% + <1 2 



(A4) 



(A5) 



Note that the well depth Vo plays no role in the development of the separable expansions as 



the product t/° 



' (o) 

Vci 



assures that it does not carry through. 



APPENDIX B: THE OPTICAL POTENTIAL 



Assuming a local form for the elastic channel element of the potential matrix, the optical 
potential for elastic scattering is defined by 

V ovt {r y. E) = Vn (r) + £ V lc (r)G { 3\r,r';E)V c , 1 (r') 

c,c'=2 

= V u (r) + AU(r,r';E) . (Bl) 
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Here AU, the dynamic polarization potential (DPP), makes this formulated interaction 
complex, nonlocal, and energy dependent as G®) are the full Green functions referring to 
the Q (= C— 1) excluded channels. Those Green functions are solutions of LS type equations 
built upon the free single channel Green function, Gq c , namely 



c 



(Q) 



(B2) 



c"=2 



This complex equation is vastly simplified when the interactions are approximated by sep- 
arable expansions of finite rank and then 0, E[ 



N 



AU(r, r'\ E) = £ Xln (r) [A(E)) nn , (E) Xn ' 1 (r') 



n,n'=l 



where 
involves 

g?\e) 



A(E) = \rj-G^\E) 



T] 



-1 



open 

E 



00 XcnWXcn'W d ° Sed ^ 



k 2 — x 2 + ie 



-x 



! dx — V] 



Xcn{X) Xcn i[X) 2 j 

h 2 + x 2 X 



(B3) 



(B4) 



(B5) 



where k c and h c are respectively. the wave numbers relevant for each open and closed channel 
considered. In application 0,0]) various perturbative- iterative schemes have been used to 
evaluate the DPP and thence to determine the (elastic) scattering S matrix. 



APPENDIX C: LOCAL POTENTIAL MATRIX ELEMENTS 

Using the Woods-Saxon forms 



1 + e^~' 



-i -l 



9{r) 



ldf(r) 



r dr 



; W ls = 2V ls Xl 



with K[ = /3 L R /a and the definitions 

T-Ro 



e(r) = exp 
A(r) 



e(r) [1 + e(r)] = — — — 



; / (r) = [l + e(r)]- 1 

= -ar g (r) 



B(r) = e(r) [e(r) - 1] [1 + e(r)] 3 
e(r) 



x 2 d 2 /(r) 




/ a 2 r\ dg 


/ de 2 





\R^) ~de 



C(r) 



[l + e(r)r 



e (r) - 4e(r) + 1 



-ar 



a \ 2 d 2 g 
Ro) de 2 



(CI) 



(C2) 



the potential matrix elements, Eq. ()37jl. take the form, omitting the relevant parity indices, 
K'c(r) = {[V + 1(1 + 1) ^] [/o(r) + 1/(8tt) «|B(r)" 
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-— w ls 

ar 



'A(r) + 1/(8tt) K 2 L C(r)] [£-s] c } 5 dc 
+V SS [f (r) + 1/(8tt) /c|S(r)] [s-I] cV 

+ + + + *(* + !)]} 

x {[4tt/(2L + 1)]* « £ A(r) [Y L -Y L ] c , c + ±(2L + l) 2 «|S(r) 



v 1 

X fe(2^+1) 



|(L0L0|£0)| 2 [Y r Y £ ] c , c 



2ar 



x [4tt/(2L+1)]5 KLJ B(r) [Y L -Y L ] c , c + 



(2L + 1) 5 



«?C(r) 



X S(2£+1) 



|<L0L0|£0)| 2 [Y r Y,] c , c 



; y ss [47r/(2L + l)]2 KL A(r) 

><E{[»-W [Y L -Y L ] C „ C + [Y l -Y l ] c , c „ [s-I] c „ c } 



1 T/ (2L + 1) 



+ 2 K ^ 



-K/ 



X 



£ {[s-I] cV , [Y r Y,] c „ c + [Y r Y £ ] c , c „ [s-I] c » c } 



(C3) 



For the specific case considered, that of neutron scattering from 12 C and allowing coupling 
with the 0+ (ground state), the 2+ (4.4389MeV), and the excited 0+ (7.64 MeV) states of 
12 C, the L = 2 multipole only is needed in the deformation expansions. Using the specific 
values for the parity Clebsch-Gordan coefficients of 



(2020|00) = -i= ; (2020|20) = ; (2020|40) = 



gives the result 

V dc {r) = {[Vo + 1(1 + 1) V u ] [Mr) + 1/(8tt) K 2 2 B(r) 



W ls 

ar 



A(r) + 1/(8tt) K 2 C(r)] [£-s] c }5 c 



+V SS [Mr) + 1/(8tt) K 2 2 B{r)] [s-I] £ 
+ {^o + ^W + !) + *(* + !)]} 



X < 

1 

2ar 



(4n/5)K 2 A(r) + -^B(r)j [Y 2 -Y 2 ] c , c + ^S(r) [Y 4 -Y 4 ] 

W ls + [e*] c } 



(C4) 
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x K/(V5)« 2 5(r) + -fC(r) [Y 2 -Y 2 ] c , c + -^C(r) [Y 4 -Y 4 ] c , c 



+ 2 V » { IW5M(r) + f B(r) J 

x£[MW [Y 2 -Y 2 ] c „ c + [Y 2 -Y 2 ] c , c „ [s-I] c » c ] 

c" 

+^S(r)X)[[s.I]^ [Y 4 -Y 4 ] c „ c + [Y 4 -Y 4 ] cV ,[s-I] c « c ]| . (C5) 

In deriving Eq. (|C5J) from Eq. (|37|l . use has been made of the matrix elements: 

(£.*) = (£'j'l'J\£ • i\£jIJ) = SwSjfSnAi + 1) , (C6) 

and 

r i if j = £ + 1 

(«•*) = x Ifti _ J- ( C7 ) 

I 2 > it J — * 2 

The spin-spin matrix element is a little more complicated |19], namely: 

(s-i) = (-)(^' + j)^f i i|(j'||i||j) (s || s || s ) 



= 5//'^(-) (1/2+J - /+/+J+ ' ) v/(2j + l)(2j' + 1)(2/ + 1) 

(C8) 

Since the operator is diagonal in I and £, and null if either I or I' is null, for our calculations 
the following particular case is useful: 

(£j'2J\ S -I\£j2J) = —L—K(l,j,f, J). (C9) 

The values of interest for K are the following: 

11 27 

K{£, £ +-,£ + -, J) = -£(£ + 2) + J( J + 1) - —, (CIO) 

1 1 23 
K(£,£- -,£--, J) = £ 2 -j(j+i) + - (Oil) 

K{£, £ + -,£ - -, J) = K(£, £--£+- j) (C12) 



We give, finally, the matrix elements of the scalar product of two rank L spherical harmonics 

as 

(Y L -Y L ) = (I'j'l'J Y L (f)-Y L (f ) IjlJ 
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(-) 



b+r+J) 



j+I'+V 



f J L 
I I' J 



l'\ ) f 



Yr(r) 



J-)j)(J* Y L (T) 



(2j + l)(2f + 1)(2/ + 1)(2/ + 1)^(2L + 1) 



{'. 


J 


5} 


/'in 




I 




1/ -t cj 



then, on using the identity 

{/ I v } {l0L ° m = H^" +l) 

Eq. (|C13|) reduces to 

'/y/'j 



(2Z + l)(2j' + 1) 



j-LO 



J 2 



(C13) 



(C14) 



Y L {f).Y L {T)\ljU) = 2^(2/ + l)(2j + l)(2j> + 1)(2L + 1) 



x(/0L()|/'())(ji/-i|L0 



i r j 



(C15) 
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